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Abstract 

The framework to describe natural phenomena at their basics being 
quantum mechanics, there exist a large number of common global 
phenomena occurring in different branches of natural sciences. One 
such global phenomenon is spontaneous quantum decay. However, 
its long time behaviour is experimentally poorly known. Here we 
show, that by combining two genuine quantum mechanical results, it 
is possible to infer on this large time behaviour, directly from data. 
Specifically, we find evidence for non-exponential behaviour of alpha 
decay of 8 Be at large times from experiments. 

PACS numbers: 23.60. +e, 03.65.Ta, 03.65.Nk 

The decay law in quantum mechanics is necessarily non-exponential [1] 
and deviates from a simple exponential form at small and large times (for 
reviews see [2, 3, 4]). Recently, some experimental facts regarding the short 
time behaviour have confirmed the expectations [5]. Till today this is not 
the case for the long time evolution. Although this long tail of the survival 
probability of an unstable quantum state is more often than not, not directly 
observable [6, 7, 8] (indeed, we are not aware of any such successful experi- 
ment), the information about resonance time evolution, in general, and the 
long time behaviour, in specific, should as a matter of principle be encoded 
in the resonant scattering data A + a — > resonance — > B + b. Since energy is 
a dual variable to time, we can imagine that the survival amplitude 

A 9 {t) = (y\e- iHt \^) = (*|*(*)) (1) 

(related to the survival probability, Py(t) = \A^ (t)| 2 /|*4*(0)| 2 ) can be writ- 
ten as a Fourier transform of an energy dependent quantity which can be 
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constructed on the premises of the scattering matrix, S. If this is so and the 
S- matrix itself for a resonant reaction is extractable from the experiment, 
the information encoded in the scattering allows us to infer on the time evo- 
lution of the resonance produced as an intermediate state in the process. To 
see this point clearly, let us first consider the method for the recasting of the 
survival amplitude as a Fourier transform of a spectral function. Ever since 
it was derived by Fock and Krylov [9] it has been the basis for most of the 
investigations on quantum unstable systems. 

Consider a resonance R* formed as an intermediate unstable state in a 
scattering process such as A + a — > R* — > A + a. Since the unstable state, 
|\&), cannot be an eigenstate of the (hermitian) Hamiltonian, we expand it 
(assuming a continuous spectrum) in terms of the energy eigenstates of the 
decay products A and a as follows: 

= [dE a(E) \E) (2) 

where \E) is the eigenstate and E the total energy of the system A + a. 
Substituting now for |\&) in (1), we get, 

A*(t) = J dE'dE a*(E')a(E)(E'\e^ iHt ^\E) (3) 
= J dE' dE a* (E') a(E) e^ tEt) 5(E - E') 

The proper normalization of \I/ tells us that \a(E)\ 2 should have dimension of 
(l/i?) and hence can be associated with density of states. In the Fock-Krylov 
method, Ay(t) then reads 

roo 

A*(t)= dEp 9 (E)e- lEt (4) 

-'-Eah. 

where E t h_. is the minimum sum of the masses of the decay products [10]. 
Since we have expanded the unstable state in terms of the eigentstates, 
the spectral function p 9 (E), is a probability density to find the eigenstates 
\E) in or, in other words, it is the continuum probability density of 
states in a resonance. Hence we can write, 

M £) = ^MS H (W. (5) 
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If we succeed, in a second step, to extract p 9 (E) from experimental data, 
we can process this data by equation (4) to conclude, now from experiment, 
on the long tail of the survival amplitude of the resonance under considera- 
tion. This would then be the closest we can ever come to pinning down the 
large times in a decay directly from experiment albeit by an indirect method. 

Theoretically, many different forms of p 9 (E) are available and many have 
been used in a rather ad-hoc fashion. The then sometimes specious results, 
especially with regard to the large times, are not always compatible with each 
other. Secondly, mostly these are purely theoretical constructs for which it 
is hard to find a bona fide connection with experiments. In all this, one 
important result originating from statistical mechanics, seems to have been 
overlooked, at least in connection with unstable states. In calculating the 
second virial coefficients B and C for the equation of states in a gas, pV = 
RT[1 + B/V + C/V 2 + ...], Beth and Uhlenbeck [11] (the derivation of their 
result is reproduced in [12], see also [13, 14]) found that the difference between 
the density of states with interaction, ni, and without, n\°\ is given by the 
derivative of the scattering phase shift 8\ as, 

n \ (°)n \ 2/ + 1 ^(E'cm) 
ni{k C M) ~ n\ '{kcM) = jr- (6) 

where / is the angular momentum of the I th partial wave and kcu and Eqm 
are the momentum and energy in the centre-of-mass system of the scattering 
particles, respectively. If a resonance is formed during the scattering process, 
Eqm becomes the energy of the resonance in its rest frame. Certainly, the 
density of states and the probability density in (5) are connected. Note 
that in the absence of interaction, since no resonance can be produced, one 
would expect (5) to be zero. Switching off the interaction (by, say, letting 
the coupling constants go to zero), n\ will not become zero, but tend to rif^ 
from above. Therefore, as long as rii — nf^ > (this is always the case for 
an isolated resonance), we can write for the continuum probability density 
of states of the decay products in a resonance, 

dProb^ffcM) _ congt d5i(E CM ) ^ 
dEcM dEcM 

which is the sought after connection between data, here in the form of Si, 
and the survival amplitude in (4). This method is a general (i.e. without any 



further restrictions with the exception of our belief in quantum mechanics) 
feasible tool for studying the time evolution of resonances from data, only 
if there are no overlapping resonances, i.e., if we handle one isolated reso- 
nance. In reality, this is difficult to realize and in most cases overlapping 
resonances will "distort" dSi/dEcyi, which can then have several maxima 
and minima, even negative (see Fig.l). These negative regions are bound 
to appear between resonances as noted by Wigner long ago in a different 
context [15]. The realistic situation of several overlapping resonances implies 
that the identification (7) is operative starting from threshold and extending 
over one resonance region, but often not beyond, as far as real experiments 
are concerned. However, one very useful feature remains when we restrict 
ourselves to large times. Large times correspond to small energies, which 
implies that in order to experimentally extract information on this region, 
all we need to know is Si at threshold and in the vicinity of the resonance. 
The exact form of how p^{E) falls off at large E, well beyond the resonance 
region, is not important to conclude on the behaviour of A^(t) as t — > oo. 
This is the reason why we restrict ourselves to large time behaviour in apply- 
ing the method of Fock and Krylov with the inclusion of the result of Beth 
and Uhlenbeck. 

To demonstrate the connection between the spectral density and phase 
shift derivative as in (7), we quote a simple example below. An amplitude 
which describes the resonant scattering process around the pole is often taken 

as, 

T = T/2 (8) 

En — E — iT/2 ' 1 ' 

from which one easily gets [16], 

dS T/2 



dE (E R - Ef + in ' 



(9) 



i.e., the Lorentzian (Breit- Wigner) form. The right hand side of (9), up to a 
constant, is commonly taken as the spectral function in (4) to display the fact 
that a one-pole approximation as in (8) leads to the exponential decay law 
[17, 18, 19]. Note that we obtained this spectral function as a derivative of 
the scattering phase shift. If a resonance, R*, is produced an as intermediate 
state in the scattering of two particles, A and a as, A + a — > R* — > A + a, 
then the energy derivative of the scattering phase shift S for this reaction, 
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gives the spectral function for the unstable state or resonance R*. If we 
compute (9) in the rest frame of the resonance, then the energy E in (9) is 
the energy available in the centre of mass of the scattering particles A and a. 
In the present work, we use the energy derivative of experimental scattering 
phase shifts to obtain the spectral function corresponding to the intermediate 
unstable state in scattering. 

Suppose now, that we have a fit to the data of 5i with a reasonable scan of 
the threshold/resonance region which will determine the large time behaviour 
and which to know is therefore a conditio sine qua non. We said already 
that in calculating the survival probability for large times, the large energy 
behaviour of every single resonance, in the situation where several resonances 
are overlapping is not of importance. We now make this statement more 
precise gaining as a byproduct more insight into the late time domain. Since 
the phase shift Si has a threshold behaviour, so will p 9 . We can parametrize 
it without loss of generality in the form p 9 oc (E C m — -E t h.) 7 to account for 
the threshold. Hence we have 

P*,(£cm) = Q{Ecu){E C M - £th.) 7(0 • (10) 

We impose the following condition on G(Ecm)' (i) G{Eth.) ^ since we 
have factorized the threshold already, (ii) G(Ecm) — > sufficiently fast as 
Eqm — > oo (in case of several overlapping resonances this is a theoretical 
assumption which, however, is inherent in the Fock-Krylov method) and (iii) 
mathematically, we allow the function to have poles zoi in the complex plane, 
i.e. l/G(zoi) = such that Qm(z 0i ) < and Ke(zoi) > 0; since the derivative 
of the phase shift carries the information about the poles [20, 21] we would 
expect on physics grounds a single pole at Er — iY/2 signifying the resonance 
parameters. This generalizes the simple Breit-Wigner form in (9) which has 
the pole but no threshold behaviour. These general properties allow us to 
compute the survival probability by going to the complex plane (though this 
method of finding the survival probability is standard and known [2, 4], we 
describe it here as it is not exactly equivalent to that in [2, 4, 22]). We choose 
the closed path Cr = C% + Csr + C 1 / 4 , starting from zero (after change of 
variables y = Eqm ~ E t h) along the real axis (Csr) attaching to it a quarter 
of a circle with radius R (C 1 / 4 ) in the clockwise direction and completing the 
path by going upward the imaginary axis up to zero (Cg). In the integral we 
let R go to infinity noting that along C 1 / 4 ^, the property (ii) and the fact 
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that e tzt oc e^ mz (Quiz < 0), ensures that the integral is zero. We subtract 
the contribution along the imaginary axis. This gives 

A^(t) = Al(t) + Al(t) (11) 

with 

A% (t) = e- iE ^ lim / dze- izt z*g(z + E th ) = C x e ~ iERt e~ r ^ 2t (12) 

; R^oo Jc R 

by Cauchy's theorem and 

Al,(t) = C 2 e- iEtht dxe- xt x>Q{-ix + E th _) 
1 Jo 

* C7 2 e-^^(^)r( 7 + l)^ (13) 

for the integral along where the approximation is valid for large times 
t. In the above, T(x) is the Euler's gamma function and C 1 is a constant 
easily calculable in terms of the parameters Er, F and E t ^ and C 2 is (— 
Equations (12) and (13) show the general features which we would expect: an 
exponential decay law followed by inverse power law corrections. The latter 
is independent of the details of Q(E CM ) as claimed, displaying also nicely 
the dual nature of time and energy. Formulae (12) and (13) substantiate 
our previous claims about the method to extract information on large time 
behaviour of resonances even if the data have several overlapping structures. 
To be mathematically precise, we note that the way we derived (13) is valid 
only for integer values of the exponent 7. We refer the reader to [4] for the 
general case, where, however, a similar result with the survival amplitude 
being proportional to l/t 7+1 is obtained for large times. 

To make a direct contact with experiment we need experimental data on 5i 
which, parametrized, can be used to extract the survival amplitude. However, 
the data have to start very close to the threshold. This is not always easy to 
find in literature, as threshold regions are not the most interesting regions to 
focus on in an experiment. 

We opted for an experiment with many data points at threshold and 
relatively small error bars. It is the a-a D-wave resonant scattering in nuclear 
physics [17-21] 

a + a -> 8 Be(2 + ) -> a + a . 
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In Fig. 1 we display the phase shift and its derivative over a wide region, 
using a simple polynomial fit to the phase shift. We find the established 8 Be 
levels, shown in the figure. Motivated by a Lorentzian form [28] with an 
energy dependent width T(E CM ), we parametrized the data in the region of 
the first 2 + resonance, in the following form: 



Si(Ecm) = tan 1 



r(E, 



CM 



En — E, 



CM 



(14) 



with 



T(-Ecm) = r 



EL 



El 



El. 



re/2 



(15) 



which is valid for the elastic case. The derivative of this parametrized 
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Figure 1: D-wave phase shifts (upper half) in a-a elastic scattering from refs 
[17-21], polynomial fit to these data (solid line) and the derivative of phase 
shift (lower half) calculated from the fit showing all established 8 Be(2 + ) 
levels, as a function of the excitation energy E ex = Ecu — -^ 8 Be(groundstate) 
and plotted here in arbitrary units (arb. units). The negative region in 
derivative of phase shift (lower half) between 5 — 15 MeV, due to the slowly 
falling phase shift is not obvious in the plot due to the scale of the vertical 
axis. 
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Figure 2: D-wave phase shift (upper half) and its derivative (lower half) in 
a-a elastic scattering as a function of Eqm ~ -^threshold, i n the region of the 
first 2 + level of 8 Be. The dashed line shows the fit mentioned in the text. 
The inset displays the accuracy of the fit near the threshold energy region 
which is crucial for the large time behaviour of the decay law. 

phase shift satisfies our requirements. In particular, k = 27 + 2 and the 
pole is highlighted by a peak structure (see Fig. 2) which is expected [20]. 
To be unbiased, we fitted k, (3, T Q and E simultaneously taking the error 
bars into consideration. Our best fit gives k = 6.36, (3 = 0.00359 GeV -1 , 
T = 0.0009 GeV and E = 7.45838 GeV and is shown also in Fig. 2. In 
the fitting procedure, special attention was paid to the mandatory threshold. 
We can take the parametrization, the derivative of the same and perform 
numerically the integration to obtain the survival amplitude. The numerical 
result for the survival probability is depicted in Fig. 3 and is our main re- 
sult. However, since we see our theoretical conditions (i)-(iii) on the spectral 
function confirmed, it suffices to use the fit together with the mathematics 
developed in formulae (12) and (13). We then conclude that at large times, 
the survival probability of the unstable 8 Be(2 + ) state at 3.04 MeV excitation 
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energy, behaves as, 

Ps Bc(t)~^. (16) 

Theoretically we would expect that near threshold [28, 29], tan<5/(i?cM) ~ 
fccy 1 , which implies, d5i(E CM ) / dE CM ~ (E CM — E th ) l ~ l l 2 also near thresh- 
old. This amounts to saying that k is expected to be 21 + 1. For the 8 Be(2 + ) 
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Figure 3: Survival probability P(t) of the decay 8 Be(2 + ) — > a + a, as a 
function of number of lifetimes after decay. P(t) = {A^^t)^ /{A^^O)^ is 
evaluated numerically using dSi/dEcM of Fig. 2 (which has been obtained 
from the fit to phase shift data) as the spectral density p 9 (E) in (4). The 
dashed line (e~ r *) shows that the decay law for the 8 Be(2 + ) state (solid line) 
is exponential up to about 30 lifetimes after which it proceeds as t -6 36 . 
and T are the resonance mass and width respectively. 

resonance, one then gets an inverse power law t~ 5 for the survival probability. 
The data on the phase shift do not seem to follow the standard threshold 
behaviour and hence we get (16). The discrepancy, however, does not look 
serious. Indeed, re-calculated the "/"-value of the fitted k is 2.68. Interest- 
ingly, the exponent 6.36 is close to the theoretical prediction of 7 for the case 
of I = 2 made in [2]. The deviation from the exponential decay law starts 
around 30 lifetimes after the onset of the decay. By this time, one could say 
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that the sample with which one started has depleted by about 13 orders of 
magnitude (~ e~ 30 ), making a direct measurement of such a phenomenon 
not feasible. The above is the case of a strong decay with short lifetime. In 
the case of weak decays, the onset of the non-exponential law at large times 
is expected to be much later [2, 3] making the direct measurement even less 
feasible. 

In summary, we combined the Fock-Krylov method to calculate the sur- 
vival amplitude of an unstable state in terms of a Fourier transform of a 
spectral function with a result in statistical mechanics of Beth-Uhlenbeck 
which identifies the continuum density of states with the energy derivative of 
the two body scattering phase shift (6) being proportional to the continuum 
probability density of states (5). Hence, using experimentally determined 
phase shifts, the method allowed us to compute the non-exponential long 
time behaviour of unstable quantum systems directly from data. An ex- 
plicit example was given and the inverse power law behaviour at large times 
determined from data. Asked as to why the merger of the two quantum me- 
chanical results, which when combined give insight into the quantum decay, 
has been overlooked so far, we can only speculate by answering that even 
the old results by Eisenbud [30], Wigner [15] and Smith [31] concerning the 
phase shift derivative have been neglected for a long time and came back into 
vogue only recently [30-34]. 

We close by quoting from [37]; "Thus it seems unlikely that nuclear de- 
cays will show deviations from the exponential decay law which they made 
famous." We have shown that it is possible, as the information on the time 
evolutions is encoded in the scattering data. 
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